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The three-body energy-dependent effective interaction given by the Bloch-Horowitz (BH) equa- 
tion is evaluated for various shell-model oscillator spaces. The results are applied to the test case 
of the three-body problem ( 3 H and 3 He), where it is shown that the interaction reproduces the 
exact binding energy, regardless of the parameterization (number of oscillator quanta or value of 
the oscillator parameter b) of the low-energy included space. We demonstrate a non-perturbative 
technique for summing the excluded-space three-body ladder diagrams, but also show that accurate 
results can be obtained perturbatively by iterating the two-body ladders. We examine the evolution 
of the effective two-body and induced three-body terms as 6 and the size of the included space A 
are varied, including the case of a single included shell, Ahu — Ohu. For typical ranges of 6, the 
induced effective three-body interaction, essential for giving the exact three-body binding, is found 
to contribute ~10% to the binding energy. 



I. INTRODUCTION 

Techniques in popular use in the nuclear three-body 
problem include the Faddeev|l|, Green's Function Monte 
Carlo (GFMC), and correlated hyper-spherical harmon- 
ics expansion methods (see, for instance, refs.0, EJ El)- 
These methods have been used with realistic phenomeno- 
logical potentials (e.g. Avl8 5j,, CD-BonnQ, etc. . .) 
with strong short-range interactions, yielding binding en- 
ergies accurate to within O.l%0. 

All of these methods treat the bare NN interactions 
in Hilbert spaces that are effectively infinite. Because 
the complexity of the Hilbert spaces grows very rapidly 
with nucleon number A, the extension of such methods to 
heavier systems becomes increasingly difficult. This mo- 
tivates another approach: Solving A-body problems in 
more tractable finite "included spaces," while accounting 
for the missing physics of the "excluded space" through 
an effective interaction. For shell-model-inspired effective 
theories that define included spaces in terms of harmonic 
oscillator (HO) bases, the missing physics includes both 
high- momentum and long- wavelength interactions. The 
former are missing because the included space contains 
only low-energy oscillator shells, while the latter are con- 
nected with the finite extent of the basis states, which 
are overconfined in an HO potential. 

Shell-model calculations often employ two-body effec- 
tive interactions, but these are generally determined phc- 
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nomenologically. In recent years there has been growing 
success in calculating these effective interactions directly 
from the underlying bare interaction in a precise and 
systematic way (e.g. refs.[3, El El EI and references 
within). The purpose of this paper is to extend upon 
this work by calculating the more complicated H e f f that 
could be used in a shell-model-inspired effective theory, 
through three-body order in the excluded space.. The in- 
clusion of three-body contributions is generally believed 
essential in building in the density dependence crucial to 
saturation. 

The exact effective interaction for the A-body problem 
is an A-body operator. The calculation of this interac- 
tion is clearly as difficult as solving the original A-body 
problem in an infinite Hilbert space. However, there are 
reasons to hope that such an A-body calculation might be 
unnecessary. This hope depends on the plausible notion 
that the effective interaction might converge as an expan- 
sion in the number of nucleons interacting in the excluded 
space. Were such an expansion to converge at the three- 
or four-body cluster level, then standard three- or four- 
body methods could be used to treat the excluded space. 
These more complicated but still tractable effective in- 
teractions could then be diagonalized in the shell-model- 
inspired effective theories to produce accurate binding 
energies and wave functions. 

There are a couple of arguments that such a scheme 
might converge. One is the success of the shell model 
when phenomenological two-body effective interactions 
are used. This suggests that a good part of the excluded- 
space physics is two-body. Another is the qualitative 
argument that short-range clustering is increasingly un- 
likely as the size of the nucleon cluster increases: such 
spacial clustering is necessary for strong, multi-particle, 
high-momentum interactions. Finally we know that the 
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Pauli principle forbids short-range s-wave interactions 
above A=4. 

Additional support for this notion comes from re- 
cent large-basis ab initio no-core shell-model (NCSM) 
calculations |j| [^. These calculations use a large 
but finite HO Hilbert space. A Lee-Suzuki (LS) simi- 
larity transformation|13j coupled with a folded-diagram 
sum|l4) is performed on the bare NN interaction to yield 
an energy-independent effective interaction. Early calcu- 
lations utilized only the effective two-body interaction as 
derived from these transformations. Hence calculations 
in large Hilbert spaces, or included spaces, were usually 
needed to ensure convergence |26j of the binding energy 
(i.e. the included space had to be large enough such that 
the contributions to the binding energy from the effec- 
tive three- and higher-body terms were small). Recent 
calculations that now include induced three-body inter- 
actions derived from the LS procedure show significant 
improvement in the convergence of binding energies for 
the alpha particle and some p-shell nuclei Il5j]. 

An alternative approach to calculating the induced 
three-body effective interaction was developed in Ref . |l(J 
using an excitation-energy- dependent three-body inter- 
action based on the self-consistent solution of the Bloch- 
Horowitz (BH) equation 16] . This method was success- 
fuly used to calculate the binding energy of 3 He. How- 
ever, both of these works to date (i.e. LS and BH proce- 
dures) involve handling the excluded-space physics in a 
very large but still finite Hilbert space. This approxima- 
tion introduces an extra dependence to the three-body 
binding on the size of the excluded space. Such depen- 
dence can only be removed by a true summation over 
all excluded high-energy modes (i.e. excluded space of 
infinte size). Though this dependence is small due to 
the fact that results of these previous works are quite 
accurate (both groups examined the convergence of their 
results as a function of the size of the excluded space), the 
rate of convergence clearly depends in detail on how sin- 
gular the underlying NN interaction is. Here we develop 
a method for calculating the effective three-body interac- 
tion that introduces no such high-momentum cutoff, but 
instead truly integrates over all high-energy modes. The 
accuracy of our results is only limited by our numerical 
precision. Thus it can be applied with equal confidence 
to any underlying NN reaction, regardless of that poten- 
tial's description of short-range physics. This interaction 
is used in the simplest test cases, 3 He/ 3 H, and the re- 
sults are compared with those that would result from 
evaluating the effective interaction only at the two-body 
level. Our treatment is based on solving the BH equa- 
tion as well, which produces a state-dependent (and thus 
excitation-energy-dependent) Hermitian interaction, 



H eff (E) = P{H + H——QH}P, (1) 



where 

A 

h = E^ + h?) 

i<j 

A A p2 

Here the intrinsic Hamiltonian H is obtained by summing 
the relative kinetic energy and potential energy operators 
Tij and Vij over all nucleon pairs i and j. The relative 
kinetic energy is found by subtracting the center-of-mass 
(CM) kinetic energy from the total kinetic energy ob- 
tained by summing Tj over all nucleons. Ma is the total 
mass of the A-body system, Q is the excluded space pro- 
jection operator, and P = 1—Q is the included-space pro- 
jection operator that defines the finite, low-energy space 
in which a direct diagonalization will be done, once H e f t 
is obtained. As E is the desired eigenvalue, which is 
not known a priori, the equations must be solved self- 
consistently. In the calculations reported here, the Avl8 
potential is used as the bare NN interaction and no bare 
three-body interaction is included. 

Because our application is to 3 He/ 3 H, rather than to 
a heavier nucleus, we will have solved the effective in- 
teraction at the A-body level. Thus we are not test- 
ing the assumption of a cluster expansion here. Instead, 
our focus is on the technique for solving the effective in- 
teraction, and its relation to the approximate two-body 
result. Using techniques familiar from Faddeev calcula- 
tions, we first determine an effective two-body-like inter- 
action, V^^f p from the scattering t-matrix. We find 
that the induced effective three-body interaction can be 
evaluated perturbatively by iterating on Vj^^ff F° r a 
certain range of oscillator parameters 6, convergence is 
achieved even for very small included spaces, including 
the limiting Ohu space. 

We also develop a non-perturbative method. By re- 
arranging the terms of the original expansion for the 
three-body interaction, the series can be summed exactly 
through numerical solution of an integral equation. This 
method is very accurate and efficient for any included 
space and any oscillator parameter. 

In Section [H] we g eneralize the two-body formulation 
presented in Ref. [Tl| . Section IIIII describes our pertur- 
bative treatment of the effective interaction, while Sec- 
tion llVl presents the nonperturbative solution via an inte- 
gral equation. The associated discussions focus on qual- 
itative issues, with derivations reserved for the appen- 
dices. We discuss the results and the need for a simple 
test of the cluster assumption (by applying current re- 
sults to 4 He) in Section Ivl 

II. PRELIMINARIES: TWO-BODY REVIEW 

To guarantee that H e ff, like the bare H, is transla- 
tionally invariant one can work in a complete basis of 
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three-nucleon HO Slater determinants that includes all 
configurations with quanta Nhto < Khuj. Such a basis, in 
traditional independent-particle shell-model coordinates, 
is overcomplete, consisting of subspaces characterized by 
the eigenvalues of the center-of-mass Hamiltonian Hqm 
that do not interact via H e ff. For A <4, it is convenient 
to avoid this overcompleteness by using a Jacobi basis, 
which reduces the A-body intrinsic-Hamiltonian problem 
to an (A-l)-body problem. This is the choice we make 
here. The Jacobi included-space (P) basis corresponding 
to the HO shell-model basis described above is thus sim- 
ple, consisting of all relative-coordinate configurations 
with NTvuj < ATiui. (We note that for A >4, the single- 
particle basis is the more efficient basis for many-body 
calculations. However, an H e ff developed in a Jacobi 
basis can be easily transformed to single-particle coordi- 
nates for use in shell- model-inspired diagonalizations.) 

Such included spaces for the two-nucleon system are 
simple to construct. For example, a A=2 included space 
consists of spatial wave functions of the form \n = 0, 1 — 
>, 1 7i = 1,/ = >, and \n = 0,1 = 2 >, coupled to 
spin and isospin wave functions to maintain antisymme- 
try. This requires / + s + t to be odd. The correspond- 
ing construction of fully anti-symmetrized A-body wave 
functions is less trivial 0|- Here we use the codes de- 
scribed in Ref. [3] to generate these for A=3. In Table |U 
dimensions of various A=3 included spaces are given as 
a function of A. 

For the three-body system one can obtain accurate 
results by diagonalizing the bare H in a Jacobi basis, 
provided A is made sufficiently large. This was done in 
Ref.[icj for the Avl8 potential with A=60, resulting in a 
binding energy accurate to ~ 20 keV. The motivation for 
solving this problem for small A using effective interac- 
tions is to explore techniques that might be more feasible 
for larger A, where the model-space growth analogous to 
Table [I] will be much steeper. 

One goal of the current work is to find techniques that 
might make the integration over the excluded space more 
tractable. In our earlier work on the two-body svstem[TH 
we found that this integration is difficult at both long 
and short distances. We found it convenient in Ref.fllT] 
to remove the long-distance difficulties, the pathologies 
associated with the overbinding of the HO, at the outset. 
As the overbinding of the HO becomes arbitrarily large as 
r increases, part of the tail of the wave function remains 
unresolved in any finite-basis treatment- though numer- 
ically the contribution of the tail becomes increasingly 
unimportant in direct diagonalizations as more shells are 
added. (This contrasts with the short-range problem, as 
Avl8 and other modern potentials are regulated at short 
distance by some assumed functional form. Such poten- 
tials can be fully resolved in finite bases, provided A is 
larger than the scale implicit in that functional form.) 
In Ref.|ll| the long-range behavior was handled by the 
following rearrangement of Eq. lfT|) 



where 



Heff = P 



E , E 

T e ff + Veff 



E-TQ 



E-QT 



P, (2) 



T eff = T + T—QT, 



Veff = V + V 



1 



E-QH 



QV. 



(3) 
(4) 



Here T and V are shorthand for J2i<j Pij anc ^ J2i<j 



respectively. Due to space restrictions Eq.[5]was stated in 
Ref. 1 1 1 1 ] without derivation. Hence we show its derivation 
for completeness in Appendix IA(27| . 

To calculate included-space matrix elements of Eq. [21 
it was first convenient to define the states 



\n > 



E 



E-QT 



P\Q >= 



E 



E-QT 



Sl>, 



(5) 



where |fi > is some state that resides in the included 
space (i.e. |f2 >£ P). With this definition, matrix ele- 
ments of Eq. |21 became 

< tt/|tf e/ /|a >= <^f\Te ff + v eff \nT>. (6) 

Our expansion of H e f t was formed by expanding the re- 
solvent of V e ff (see Eq. QJ in powers of QV: 



1 



-rQ = 



1 



-Q 



l 



QV- 



1 



-Q 



E-QH" E-QT" E-QT" E-QT 

(7) 

An iterative procedure was then used to calculate 
included-space matrix elements of this expansion. At 
each order we solved for the self-consistent energy E. 

The rearrangement of Eq. was applied to the 
deuteron in Ref. an d led to excellent results in pertur- 
bation theory for suitably chosen oscillator parameters b 
(only few orders of QV were needed to reproduce the 
deuteron binding energy to within one keV). Figure ^ 
shows that e ^q T acting on the HO state |fi > drasti- 
cally changes that state's large-r behavior: The asymp- 
totic behavior of |0 > is e _7r Ayr, where 7 = VmE, m 
being the mass of the particle As this is the proper 
fall-off, the overconfining effects of the HO are thus re- 
paired. The effects are important for the deuteron, which 
has an extended wave function because of its small bind- 
ing energy. The operator e E q T modifies only those HO 
states that reside in the last shell of the included space. 
For all other states within the included space, the opera- 
tor is identical to unity. The "endshell-corrected" states 
thus control all of the proper asymptotic behavior. 

The resummation of the kinetic energy operator al- 
lows one to adjust basis wave functions to absorb much 
of the short-range behavior of the NN interaction into 
the included space, leaving a weak residual interaction 
that can be handled perturbatively. By adjusting the os- 
cillator parameter b to low values (~ 0.4 — 0.5 fm), the 
residual Q-space contributions to the binding energy due 
to V e _}q H QV can be evaluated to within ~ one keV in 
third-order perturbation theory, even for very small in- 
cluded spaces. This is only possible because the kinetic 
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energy has been summed to all orders independent of 6, 
thus guaranteeing that the choice of a small b will not 
alter the wave function at large r. 

The optimal b presumably provides the best resolution 
of the hard core without substantially altering the ability 
of the basis to reconstruct the intermediate-range poten- 
tial. 

In the next section we will explore a similar expan- 
sion for 3 H/ 3 He. We find modifications are necessary to 
account for the disparate length scales that come about 
from the inclusion of a third particle. Before starting this 
discussion, we end this section by presenting compact ex- 
pressions for the operators e E q T P and e 1q T Q- Such 
operators appear frequently in subsequent formulae. 



A. e E qt P ano - e-qt Q operators 



The operator 



closely resembles the free parti- 



E-QT 

cle propagator. However, due to the projection opera- 
tor Q in the denominator, calculating matrix elements of 
e E q T P requires a little ingenuity. Here we show that 
the excluded-space Green's function can be expanded in 
terms of the much simpler full Green's function and op- 
erators within the included space that can be inverted 
easily. Following Ref.0]> we first consider 



PT- 



E 



E-QT 



P = PT{ 



E 



-P- _ 1 _PT _ } P, 



E-T E-T E-QT' 



(8) 



where Q = 1 — P. Collecting terms in this "Schwinger- 
Dyson" form and noting that P 2 — P gives 



P<l + T 
=>PT 



1 



E-T 
E 



P PT- 



E 



E-QT 



-P = PT- 



E 



E-T 



:P 



E-QT 



-P 



P^P 

E-T 



PT - P. 

E-T 



The last term in Eq. |UJ PT E E T P, can be rewritten 
PT - P = EP{ P - P -1)P. 



(9) 



E-T \ E-T 

Plugging the above equation into Eq. [5] finally gives 

E ( f 1 



PT- 



E-QT 



P = P { E- IP 



E-T 



-P 



P. (10) 



Similarly, the "Schwinger-Dyson" form of e _q T P is 



E 



-P 



E 



-P 



E-QT' E-T' E-T 
Using Eq. (fTTTf) then yields 



PT- 



E 



E-QT 



E 



E-QT 



-P = 



1 



E-T 



P- 



E-T 



-P 



P. 



(11) 



(12) 



This expression is relatively simple to use. The free par- 
ticle propagator, E _ T , is known analytically. Indeed, 
its form is diagonal in momentum space. The operator 
P E _ T P represents a matrix composed of included-space 
overlaps of the free particle propagator. Since we work 
within small included spaces, inverting this matrix is not 
difficult. The included-space matrix elements are easy 



to calculate, as analytic expressions for < f2/| 



E-T 



fii > 



for any A-body system exists. These expressions involve 
multiple sums over hypergeometric functions and gamma 
functions. In Appendix[U]we give the relevant expression 
for the three-body system. 

To derive an analogous expression for e ^q T Q, we first 
consider 

1 1 1 - 1 



E-QT E-T 
1 1 



E-T E-T 



E-QT' 
1 



-PT- 



P- 



E-T 



-P 



E-T 
i rp x 



(13) 



EE -T' 



where we have substituted Eq. ^| in the second line 
above. With the use of Eas.lT2land [l31 e 1q T Q becomes 



E-QT" E-QT 



1 



1 



E-T E-T E-T 



P- 



-P 



(14) 



E-T 



The physical interpretation of Eq. 1141 is simple: the term 
on the LHS represents free propagation in the excluded 
(Q) space, while the first term on the RHS (bottom line) 
represents free propagation in the full (P + Q) space and 
the second term subtracts off the contribution coming 
from included-space (P) propagation. 

It is convenient to make the following definitions: 

Go = TT~~7^ i (15) 



E-T' 



p- 



l 



-p 



1 * E-T' 

With these definitions, Ea s . 1 1 21 and 1 1 41 become 



(16) 



E 



E-QT 
1 



P = G Tv 



E-QT Q - G °- G ° r o G °- 

Finally, with the expressions above, we can express our 
perturbative expansion in a concise form, 



where 



H eff = E - r o 1 + r o lr i r o 1 

+r - 1 {r 2 -r 1 r- 1 r 1 }r- 1 - 

T x =PG VG P 
T 2 =PG a VG VG„P 

T n =PG Q (VG ) n P 



(17) 



(18) 
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In the expansion given by Eq. El E — T 1 corresponds 



to 



E-TQ ±e ff E-QT ' 



while r-^if- 1 



tV-. 



and so on. 



q corresponds to 

As we stressed in Ref . [TT| . 



s-tq ' E-QT 

matrix elements of the operators shown in Eg. 1181 can be 
simply calculated using a recursive procedure. 



III. A PERTURB ATIVE THREE-BODY 
CALCULATION 



Similarly, one might apply the expansion given by 
Eqs. (TJJ and (|T%|) directly to 3 He/ 3 H. As was done in the 
two-body system, we construct endshell-corrected states 
to build in the correct aymptotic forms. Figure [21 shows 
momentum-space results for one lOtko endshell state. As 
expected, the modified wavefunction acquires a strong 
peak near zero-momentum, corresponding to large-r cor- 
rections. 

To determine an optimal b for subsequent calculations, 
H is minimized within the included space. The results 
are given in Fig. [3] for several included spaces. As ex- 
pected, the endshell states improve the binding of the 
three-body system and lower the optimal b, relative to 
uncorrected HO results. However, the improved bind- 
ing energies are still far away from the exact answers [28). 
This contrasts with the deuterium results (see Fig. 2 of 
Ref. ^l]]) , where Oth-order results using endshell states 
were accurate ( ~ 50 keV), and became nearly exact in 
second- or third-order perturbation theory. The 0th or- 
der A=3 results underbind by ~ 3 MeV. 

As Fig. 0| shows, these results are not corrected in 
perturbation theory. Successive orders produce wildly 
oscillating values for the binding energy, in contrast to 
the rapid convergence found for the deuteron. The dif- 
ficulty is that a single distance scale b does not provide 
sufficient freedom to describe short-range behavior gov- 
erned by two relative cordinates. Extended states such 
as those of Fig. [S] are problematic. Once one makes a 
choice of Jacobi coordinates, the short-range behavior of 
some two-body cluster will be difficult to describe. For 
example, the choice (r 3 — r 2 )/-\/2 and (2fi — r 3 — r^)/v6 
prevents one from adjusting b to account for the short- 
range interactions of nucleons 1 and 2. Consequently, 
important hard-core interactions lie outside the included 
space, leading to nonperturbative behavior. 



A. Invoking Faddeev Decompositions 

Because the deuteron was easily treated, it is possi- 
ble to sum to all orders the repeated excluded space 
scatterings of any two-nucleon cluster. Such a partial 
summation should account for the strong repulsive inter- 
action between any two coupled nucleons, leaving only 
intermediate- and long-range interactions. Such interac- 
tions can be treated reasonably accurately within the in- 
cluded space using our endshell-corrected states, leaving 



perturbative corrections. This partial summation corre- 
sponds to the Faddeev decomposition 0, |2(j. The start- 
ing point is again Eq. 

< n f \H eff \m >= <%|T e// + v eff \su>, 

where we cast V e ff in its integral form, 

Veff = Veff = V + V-±—QV? f} . (19) 



E-QT 



Here the superscript 3 denotes that this is a three-body 
effective interaction. Now 

= | >, which satisfies 



effective interaction. Now consider the state V^A Q a > 



> = v\ n a > + V 



1 



E-QT 



Q\ >■ 



(20) 



This can be decomposed into Faddeev components 
>ij, given by 



>23 =^23 1 Or > + V23 

> 3 i =Vsi\ n a > + V31 



1 



E-QT 
1 

E-QT 
1 

E-QT 



Q\ >, 

Q\ *7>, (21) 

Q\ *7>, 



where 



> = I > 



12 



23 



31' 



(22) 



Equations [23 and [21 show that the Faddeev components 
require the solution of three coupled integral equations. 
For example, | ^ a > 12 satisfies 

I >12 = SI* > + 

Vl2 F^~QT Q *^ >12 + 1 *^ >23 + 1 ^ >31 

.J 23 ) 

Eauationl23lcan be inverted with respect to the | ^ a > 12 
Faddeev component, giving 



**>12 = ^12,3/1 fi* >" 



V^fl^QfQ (l 9 ° >23 + I >3l) . 



(24) 



where 



V, 



(2+1) 



12,e// - 1 _y 1 ( 
1 V12 E-QT K 



'-V 12 + Vv. 



-Vx2 
1 



E-QT- QV 12 



QV r . 



(25) 



= Vi 2 + V v . 



1 



E-QT 



QV {2+1) 



Analogous equations can be found for the remaining Fad- 
deev components. The superscript ( 2+1 ) indicates that 
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the effective interaction given by Eq. [33] represents the 
sum of repeated potential scatterings between the same 
two nucleons (in this case, nucleons 1 and 2), while the 
third nucleon remains a spectator, as illustrated in Fig. El 
This expression is exactly the partial sum mentioned 
above. Notice that the middle expression in Eq. [23] is 
very similar to Eq. Furthermore, the integral form of 
U2 ~eff ^ s sminar to that of the G-matrix[2l] . Yet there 
are subtleties that differentiate the two. We will return 
to this point later. The integral equation given by Eg. 1251 
can also be solved exactly in terms of the (2+l)-body 
t- matrices 20] . In Appendix we show 



1/(2+1) 



tl2Go 



r + rt +1 



Goti 



where 



t\2 — V\2 + Vi2Go<12 



(26) 



(27) 
(28) 



and To is given by Eq. 1161 The physical interpretation of 
Eq. |2S1 is similar to the one given below Eq. the first 
term on the RHS represents the effective interaction com- 
ing from repeated scatterings in the full (P + Q) space, 
while the second term subtracts off the contribution from 
the included (P) space. 

The advantage in using the Faddeev decompositions 
comes from exploiting particle exchange symmetries. As 
the states | U, a > are fully anti-symmetric, the individ- 
ual Faddeev components can be related to each other by 
simple permutation operators |20j. In particular, 



* CT > 23 =Pi 2 P 23 | >12, 
* CT > 31 =Pi 2 Pi 3 | *a > 12 , 

• | *T> =(i + n)| *7> 12 , 



(29) 



where the permutation operator II is defined as II = 
P12P13 + Pi 2 P 2 3- The last expression in Eq. [23 then 
becomes 



which can be formally inverted with respect to | ^ a > 12 
to give 



W " >12 — \ V 12,eff 



v, 



(3+0) 



12,e// 



to a >. 



(30) 



where 



(3+0) 



12,e// 



(2+1) 



u,eff E - QT - QIIV]^,!// 



^nvg. (3i) 



EauationlSTldefines the induced three-body interaction 
coming from repeated scatterings of V^tf} ( see Fig.[7J|. 
The operator II ensures that every insertion of V^Xf\ 



represents scatterings coming from different pairs of nu- 
cleons (i.e. the rungs of Fig. [7| alternate), preventing any 
double couting. As there are no spectator nucleons in 
this expression, it is labeled with the superscript ( 3 +°). 
Matrix elements < H e f f > are given by 



< a t \H eff \ 0/ > = 

I < m 

= <M t \T eff +3 



ff 



V, 



(3+0) 



lm,eff 



n f > 



(32) 



r, 



(3+0) \ 



12,eff) I "/ 



fit >, 



where the factor of 3 in the last expression comes from 
the anti-symmetry of the states. An expansion of Ea.1321 
can now be found by expanding the resolvent of Eq. 1311 
m powers 01 11Vy 2 e // > 

E - QT - QWjg$ ~ 

E-QT Q + E - QT QnV ^f E - QT Q + ' ' ' ' 

Expressions analogous to those of Eq. I|18|l can be con- 
structed from II and V^~lff, 



P 



r„ -PG (v^jfGo^ 

7 „ =PG (nvg$ Go) * P (33) 

f n =PG V r / 2 2 + / 1 }G {llV^f f G^ P. 

Combining Eq.[2Slwith the expansion of Eq. 1321 gives the 
final result (compare with Eq. I17fl. 

H eff = E ~ T o 1 + 3r o r iT^ + 



3ro 1 {f 3 -r 1 rj 1 7 1 }ro 1 + 
3To l {f 3 - f 2 ro x 7l - r 1 r - 1 72 + 

rir SiFq S1K0 1 + 



(34) 



Figure [HI shows the convergence of the binding energies 
for 3 H and 3 He for several included-space and oscillator 
parameter choices. The expansion of Eq. I|34|l converges 
even for very small included spaces, including A = 
(which corresponds to a single matrix element). Surpris- 
ingly, the most rapid convergence occurs for the smallest 
included spaces. We do not have a convincing physical 
explanation for this result. 

Note that the zeroth order results (i.e. E — T^ 1 + 
Sr^r^ 1 ) all overbind in Fig. 03 This is possible since 
the zeroth order calculation is no longer variational, as 
Xff now depends on E (sec Eq. [3SJ . Hence there is no 
prescribed method for finding the optimal 6, other than 
by trial and error. (We stress that fully converged re- 
sults will be independent of 6, as we have executed the 
effective theory faithfully. By an optimal b we mean one 
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that will speed the convergence.) The values of b used in 
Fig.|S]should be near the optimal. Note that these values 
are much larger than those found for the deuteron cal- 
culations of Ref.0|. Since v/H// i s calculated exactly, 
all short range two-body correlations are correctly taken 
into account. Hence b is no longer forced to small values 
by the demands of short-r physics, but instead can re- 
lax to values characteristic of the size of the three-body 
system. 

As mentioned earlier, V 12 ^| is reminiscent of the G- 
matrix found in traditional nuclear many-body theory, 
which is also an infinite ladder sum in particle-particle 
propagation. However the differences are substantial. 
Traditional shell-model calculations have always used G- 
matrix interactions at the two-body level, ignoring any 
dependence the operator may have on spectator nuclc- 
ons (including in general certain violations of the Pauli 
principle). The operator V^^ji, on the other hand, de- 
pends explicitly on the spectator nucleon kinetic energy, 
as well as on the kinetic energies of the interacting nu- 
cleons. This dependence is manifest in Eq. [23 as the 
operator T in the resolvent e ^q T Q represents the sum 
of the kinetic energies of all nucleons. This dependence 
on the spectator nucleon is essential for calculating the 
correct two-body correlations within the three-body sys- 
tem, as Q is defined by the quanta carried by the three- 
body configurations. Similarly, for A-body calculations, 
the analogous operator, V^~^f^~ 2 ^\ would carry the de- 
pendence of all spectator nucleons. 

It appears that a perturbative expansion for the bind- 
ing energy of 3 He/ 3 H results only if one first sums two- 
body ladder interactions to all orders (i.e. V^^f/)' This 
summation softens the two-body interaction, which on 
iteration can then generate the effective three-body inter- 
action (i.e. V^t//) perturbatively. A similar procedure 
was followed in Refs.j^M where the triton binding 



energy was calculated by perturbing in di-baryon fields. 
The di-baryon fields themselves represent an infinite scat- 
tering of two-body interactions, in analogy with V^tf}- 

IV. A NON-PERTURBATIVE THREE-BODY 
CALCULATION 

A drawback of the perturbative approach of the previ- 
ous section is that for each included space defined by A, 
there is only a limited range of bs for which the expansion 
converges readily. Furthermore, there is no definite pro- 
cedure for estimating the optimal b, as the zeroth-order 
calculation is no longer variational. Thus a nonpertur- 
bative procedure would be attractive: any b could be 
chosen, and physical observables would prove to be inde- 
pendent of 6, as they must for a rigorously executed ef- 
fective theory. Such a non-perturbative summation is in- 
deed possible by reshuffling the terms of Eq. (|34[1 to form 
a summable geometric series and by numerically solving 
an integral equation. This approach provides an oppor- 
tunity to directly compare the relative sizes of V^^tt 

and V^tf} as functions of both A and b. 

In this section we only show results for the triton sys- 
tem, as 3 He calculations produce similar results 29] . The 
starting point this time is Eq. 1341 

ff e// = £-r - 1 + 3r - 1 r 1 r - 1 + 
3r - 1 {f 2 -r 1 r - 1 7l }r - 1 + 

3r 1 {f 3 - f 2 ro X 71 - r 1 r - 1 72 + r^ 1 ^ 1 ^} r„ 1 
+ ..., 

(35) 

where the set {r„,r„, 7 „} is given by Ea. 1331 This ex- 
pansion can be rearranged into 



H eff - E - r o 1 + 3r o lr i r o 1 

+3(r - 1 (f a + f 3 + . . .) IV 1 - IV 1 (f a + f 3 + . . .) r„ 1 ( 7l + 72 + . . .) 1^ 1 + ...) (36) 

-3(17 ^r^ 1 (71 + 72 + ■ • •) iv 1 - r^r^ 1 ( 7l + 72 + . . .) 17 1 ( 7l + 72 + . . .) 17 1 + ...). 

1 

Equation 1361 can be verified by expanding its terms out 
and directly comparing with Eq. 1351 The following defi- 
nitions can be made to simplify Eq. 1361 

f oo = f a + f 3 + • • • (37) H *ff = E ~ r o 1 + srVriiV 

7 oo = 7 i + 7 2 + • ■ ■ , + 3 ( r o : r 00I7 1 - r^f ooTq SodY 1 + . . .) 

-3(17 x r x ro 1 - t^t.t^t^t^ + ...) 

where the above expressions represent infinite summa- (38) 
tions. Substituting these expressions into Eq. 1361 gives The two geometric series in the equation above can be 
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summed to give the final desired result, 

H eff = E - + 3F lr i F 1 



-3IV 1 foo-riroSo 



1 + To So 



o (39) 



-T _i_ QT/( 2 +!) , ot/(3+0) 



where we have made the following identifications, 



T eff — E - r 



t/(2+i) _r-ir r- 1 

V 12, eff — 1 o i o 



eff 

t/( 3 +0) j-i— x i -p 

V 12,eff 1 1 " 



'°° / 1 I T>-1 



The usefulness of Eq. (|39[) depends on having an ef- 
ficient method for calculating the included-space matrix 
elements of and r^. This can be done by first consid- 
ering the state 7oo| Q c >= >, which, using Ea. 13*71 
satisfies 

>= 7 i| n a > + 72 | n a >+... 

= G HV 12 Go\ n a > +G UV 12 G UV 12 G \ Q a >+... 
= G HV 12 Go\ fl a > +G UV 12 \^ ff > . 

(40) 

Hence \*f?„ > satisfies the Fredholm integral equation 
of the second kind 24]. There are numerous numeri- 
cal methods available for solving this particular integral 
equation. We use a general projection algorithm pi]. Due 
to the large number of partial waves and the complicated 
structure of the operator II, solving the integral equation 
by matrix inversion is impractical. 

Once l^o- > is found, the matrix elements 

< ft 4 |7oo| n a >=< > 

can be evaluated. Then the matrix elements < 
n^Pooin^ > are easily calculated since 

>=G y 12 |* ff >, 

which can be verified using Eas. 1331 and 1371 Hence 

< Oiifooi n a >=< n t \G v l2 \y a > . 



Tables iHl and Hill show the calculated binding energies 
for the triton system as a function of the included-space 
size A and oscillator parameter b. In Table [H] we have 
ignored the small isospin- violating parts of the Avl8 po- 
tential. Hence isospin is conserved at T = 1/2. Table InTl 
includes isospin- violating contributions. This causes a 
small admixture of T = 3/2 components into the ground 
state. At about the keV level, the results are independent 
of the choice of the included space, i.e., of b and A. The 
~ few keV variations illustrate the level of our numerical 
precision. Note that for the values 6=.83 and 1.17 fm, 
our calculations of the previous section would give non- 
converging binding energies. This of course emphasizes 
the non-perturbative nature of these calculations. 



Note that Eq. has the correct limiting behavior as 
A — > 00, H e ff — ► H. In this limit the projection operator 
P 1, giving T^ 1 



Gq 1 . Thus 



E-T 



— I A — >oc 



T. 



It is also straightforward to show that 



r, 



(2+1) _ 



12,e// 



ro'riro 



1 A- 



1 



1 + ti 2 Go 



12 



where we have used Eq. [23 to obtain the second line 

vanishes in this 



above. Finally, it is obvious that V^tf} 
limit since 



To 



I ll 7oo - 



G V lo 



0. 



Hence H e f *—* H. 

Tables IIVI IVII show the variation of the matrix el- 



ements of 



f and V^2~^ff with A and b. 



12,e/f a ^ r 12 eff w — " »■ ISGSpin- 

violating effects have also been ignored in these results. 
As the number of included-space examples is small, it 
is difficult to extract any limiting behavior. Yet it does 
seem that as A increases, matrix elements of V^,}\ a P~ 
proach the bare V\ 2 (indicated by the symbol 00 in the 
tables). This is especially evident in Table [V] where the 
effect of renormalization is small (i.e. the renormalized 
matrix elements are not too different from their bare ma- 
trix elements). In the cases involving V/jfe//' limiting 
behavior is also difficult to extract. Indeed, in Tables. llVl 
and IVII matrix elements of V^^ft tend to grow away 
from zero with increasing A rather than tend toward zero. 
However, in these cases the renormalization is strong, so 
that limiting behavior would not be expected for such 
small A. In Table Ivl the renormalization is much weaker, 
and the limiting behavior ZV^Xj^f ~> is clearly seen for 
the Tiuj = 15 and = 60 MeV examples. It is also clear 
from Tables. ITvT IVII that the matrix elements of V, 



(2+1) 
12,e// 

and e// have a strong dependence on A and b (Tito), a 
dependence often ignored in traditional shell-model cal- 
culations. 

Generally the contribution of < ^// > to the overall 
binding energy is small compared to < > (~ 10% 

of ~l~f \ )j a s is evident from Fig. El This graph gives the 



< T eff >, 3 < V^f'f >, and 3 < V^ft > contribu- 
tions to the binding energy as functions of A and b. Inter- 
estingly, for each included-space A, there are two values 



of b at which the contribution of < V, 



(3+0) 



12,e// 



> completely 

vanishes. Table IVIII lists specific values. If this result 
proves to be a generic property of the three-body H e ff 
in heavier nuclei (that is, if zeros appear in those calcu- 
lations at approximately the same value of b) , this could 
greatly simplify the more complicated included-space di- 
agonalizations required in heavy systems. For example, 
the three-body contribution could be ignored or could be 
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explored in low-order perturbation theory. This would 
require establishing that the A=3 zero for < > 
does persist at about the same b in heavier systems. This 
seems plausible, as the three-body ladder in a heavy nu- 
cleus is quite similar to that in A=3. 

V. CONCLUSION 

We found in Section II I II that the simple perturbative 
scheme found in Ref.[TI| (ie. Eq. $171 ) did not directly 
extend to the three-body system due to the additional 
length scale introduced by the third nucleoli. For any 
choice of Jacobi coordinates, this second length scale ex- 
cludes a class of hard-core interactions from the included 
space regardless of the choice of b. This problem was 
circumvented by invoking the Faddeev decomposition on 
H e f f , which allowed us to sum to all orders the repeated 
potential scattering between the same two nucleons, gen- 
erating Vi2~^jf- Hence the strong repulsive two-nucleon 
force was treated non-perturbatively. The remaining con- 
tribution to H e f f , Vi2*fl , can then be calculated by sum- 
ming repeated insertions of ~ljf on alternate pairs of 
nucleons. For certain ranges of o, this expansion con- 
verged after several orders of perturbation even for the 
smallest allowed included spaces, as shown in Fig. El 

In Section Hvl we described a non-perturbative method 
for calculating -ff e // in the three-body system, summa- 
rized in Eq. (|39|l . This allowed a direct comparison of the 

contributions of < T eff >, < V^j} >, and < V^f f >. 
Because the calculation was non-perturbative, we were 
able to explore the dependence of the included-space 



matrix elements (and establish the independence of the 
binding energy) for a wider range of 6s and As. Also, 
for each included space, we found two values of b where 
the total binding was due to vj^ *H , a result that should 
be explored in heavier systems. The numerical methods 
included an efficient algorithm for solving integral equa- 
tions (i.e. Eq. I|40|l). In the current applications to A=3 
systems, the numerical effort was not substantial. 

Work is in progress to apply a similar formalism to 
the alpha particle. As in the three-body case, one 
must account for the various length scales inherent to 
the four-body system. This can be done by invok- 
ing the Faddeev- Yakubovsky decompositions on H e ff, 
in direct analogy to the three-body decompositions pre- 
sented in this paper. Such a procedure will not only 
yield < T eff >, < V$$) > and < > (di- 

rect analogs of < V^^jt > and < K^f) >), but also 

< Vi2^// >• This last expression represents the induced 
effective four-body interaction. It will be interesting to 
see if these calculations verify the commonly implicit as- 
sumption of the hierarchical sizes of these interactions 
(i.e. Vg+? f < Vg+}} < Vg+f { ). The four-body 
effective interaction is the most complicated lowest-order 
(s-wave) operator that can be constructed. Furthermore, 
alpha-particle clustering is an important phenomenolog- 
ical feature of nuclear structure, apparent even in simple 
systems like 8 Be. Thus it will be interesting to deter- 
mine the size of this contribution. The A=4 calculations 
will also check whether zeros in the three-body effective 
interaction persist. If they do, we may learn something 
about their "trajectories" in b as A is inceased. 



APPENDIX A: DERIVATION OF EQ. [T7| 



For completeness, we show the derivation of Eq. [21 which was first used in Ref. 

El 

in deriving a perturbative 
expansion for the deuteron. We begin by explicitly expressing the BH equation in its various components, 



H, 



eff 



--H + H- 



1 



E-QH 



QH 



=H + T ^ -^z QT + V - -p-r- QV + T — —- QV + V ' 



E-QH 



E-QH 



E-QH 



E-QH 



QT. 



Next consider the operator 
T 1 



E-QH 



QT =T 

=T 
=T 



1 



1 



QV- 



1 



E-QT E-QT^ E-QH 



QT 



1 



=TQ 



E-QT 
1 

E-QT 
1 



l + QV 
l + QV 



1 



E-QH 
1 



QT 



1 



QV- 



1 



E-QT E-QH^ E-QT 



QT 



E-TQ 



E-QT + QV + QV 



1 

E-QH 



QV 



1 

E-QT 



QT, 



(Al) 



(A2) 
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where the relation e 1q T QT = TQ e 1 T q was used. A similar calculation gives 

1 



' QV=T ' 



E-QH 
1 



QT = 



E-QT 



V + V 



V + QV 
1 



E-QH 



QV 



E-QH 
1 



E-QH 

Substituting Ecis. IA2I and IA3I into Eci. lAll and using the fact that 

1 E 



E-QT 



QV 
QT. 



1 



E-QT 



QT - 



E - QT' 



and 



T + TQ ^ 1 (E-T) „ 1 QT = TQ ^— — (T - T^-T)- ' 



E-TQ' 



E-QT 



E-TQ 



E 'E-QT 



QT, 



gives the desired result: 



H eff = T Ql^-7F^(T - T %T)l^-7^Q T + T QlF-^( V + V v n^ W- ' 



E-TQ 



E 'E-QT 



E-TQ' 



E-QH" 'E-QT 



QT. 



Hence 



where 



< n' I' a\H e ff\n I a >= < n' I' ct\T e ff + V e ff\n I a >, 



Teff = T + T^T 



Vef J = V + V 



1 



E-QH 



QV, 



and 



I a > = 



E 



E-QT' 



I a > 



(A3) 



(A4) 

(A5) 
(A6) 

(A7) 



Note that the resolvent in Ea. lA6l is now sandwiched between Vs only, and not Ts. Forming a perturbative expansion 
involves expanding this resolvent as was done in Eq. [7] 

The operator e _^q T acts non-trivially only on endshell states. That is, 



E i i ^_ f \n I a > if \n I a > endshell 
E - QT ' \ \n I a > if \n I a > below endshell, 



(A8) 



which can be easily verified by expanding the operator in powers of QT and noting the fact that the kinetic energy 

operator T can change the principal quantum number n by at most ±1. These 'tilde' states \n I a > represent the 
culmination of the infinitely repeated scatterings of the operator QT on the original states \n I a >. Such scatterings 
persist at large-r even though the potential V has already vanished. 



APPENDIX B: DERIVATION OF V ( ntll (EQ. 1251) 



12, eff 



Iterating the last expression of Eq. [23 on vf' ' 



12,e// S lves 



Vv. 



1 



-QV 



QV 12 - 
1 



(Bl) 



E-QT" E-QT 



QV r . 
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Replacing e ^q T Q by the relevant expression below Eq. ^1 Eg. IBll becomes 

O =Vi2 + Va (Go - GoT-'Go) V 12 

+V 12 (Go - GoT^Go) Vx2 (G - GoTo X G ) V12 ( B2 ) 
+ ... 

The various terms of Eq. IB 21 can be re-ordered to give 

V 12 + e}) = (^12 + V ^G V 12 + V 12 G V 12 G V 12 + . . .) 

- (Via + V 12 G V 12 + . . .) Gor^Go (V 12 + V 12 G V 12 + . . .) 
+ (V12 + V 12 G V 12 . . .) GqT^Gq (V12 + V 12 G V 12 . . .) Gor^Go (V 12 + yi 2 G V 12 ...)-..., (B3) 

where 

V12 + Vi 2 G Vi 2 + Vi 2 G Vi 2 G Vi 2 + . . . (B4) 



represents an infinite sum of two-body interactions between particles 1 and 2. The solution to Eq. IB4I is given by the 
operator t\2, where 

*12 = V 12 + ^ 12 G ti 2 . (B5) 

As mentioned near the end of sect. IIIII the operator t\ 2 differs from the usual two-body ^-matrix since it is imbedded 
within a three-body space. Hence Go contains a dependence on the Jacobi momentum q of the third 'spectator' 
nucleon. This dependence carries over to the operator t\ 2 . Hence the designation (2 + l)-body i-matrix. 
Replacing the infinite sum of Eq. IB4I by t\ 2 and inserting this into Eq. IB3I gives 



^12 ~eff — ^i2Gor o 1 Go<i 2 

+ £i2G r Q 1 r^ +1 ^r c 

, r r -l r (2+l) r -l r (2+l) r -l r r . 



^i2G r Q 1 r£ 2 3 +1 )r 1 G ti 2 mg-j 



where r^" 1 " 1 - 1 is given by Eq. i.e. 

r£ +1 ) = PG ti 2 G P. 
Equation IB6I can now be summed geometrically giving the desired result 

Vl2,l// = ^12 — ti 2 G() To + r^, +1 ^ Go^l2- 



APPENDIX C: MATRIX ELEMENTS OF ^ FOR THE THREE-BODY PROBLEM 



The relevant matrix element is 



p 2 p' 2 



< n'l'-N'L'\——\nl-NL >= 5 Vtl 8 L ., L / dpdp 1 R* n , ^{p)R* N , >L {p')- ^ ^— R n ,i(p)RN,L(p')- (CI) 

h — 1 J h — p z / 2m — p z / 2m 

Replacing the radial integrals by their series expansion and changing to dimensionless variables, Eo. ICll becomes 



5i'jS L ,, L (-l) n+n ' +N+N ' v/Gi(n,n',/;,iV,7V',L)|- 



n,n 



N ' N ' (_i\m+m'+M+M' r , „ „2(l+l+m+m')„r2n+L+M+M') 

\^ \ tl / dxdx'e-^ +x - (C2) 

^G 2 (n, < I; N, N', L; to, to', M, W) J ff-^-x' 2 



m,m 
M } M'=0 



12 



where 



din, n', I; N, N', L) = ATin)Tin')T{N)T{N')T{n + 1 + l/2)r(n' + I + l/2)r(iV + L + l/2)r(iV' + L + 1/2), 



and 



C 2 in, n',l;N,N',L; m, m! , M, M') = 

T{n - m)r(n' - m')T(N - M)T(N' - M')T(l +m + 3/2)T(l + m' + 3/2) 

T(L +M + 3/2)r(Z' + M' + 3/2)r(m + l)T{m' + 1)T{M + l)T(M' + 1). (C3) 



An analytic solution to the integral in Eg. IC2l can be found by first changing variables to cylindrical coordinates: 

X =pCOS(j) 

x =psin4> 

2 2 i 12 
p =X + X . 



Hence the integral becomes 

2(l+l+m+m') I2(l+L+M+M') 



S 

JO 



dxdx'e- {x2+x ' 2) - 



^-x 2 - x' 2 



dp e 

7T/2 



P 2 P 



2(2+l+m+m' +L+M+M') + 1 



dip icoscf)) 2 



2(l+(+m+m') / S ^^2(l+X+M+M') 



(C4) 



The integrals in square brackets have analytic solutions that can be found in any book of integrals (e.g. see Ref . |25| ) . 
They are written here for completeness: 



dp e 



2(2+l+m+m' +L+M+M') + 1 



I^Yi+l+m+m'+L+M+M' + l f 2E\ 



2+l+m+m' +L+M+M' 



TlLO I 



2E S 



x T(3 + l + m + m' + L + M + M')T(-(3 + l + m + m' + L + M + M'), -—) (C5) 

Tito 



(i 



" /2 d^ ( CO #) 2 ( 1 +^-+-') (sin0) 2(i +i+ M + M' ) = Td + rn + m' + mnL + M + M' + m^ 

Y v Y> y YJ 2r(3 + l + m + m' + L + M + M') y ' 



This procedure generalizes for matrix elements of E ^ T for higher-body systems. For an A-body problem within 
Jacobi coordinates, a change of variables to (A-l)-dimensional spherical coordinates is needed so that the resulting 
(A-l) integrals 'factorize' (e.g. those in square brackets of E q. IC4|I. 
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FIG. 1: Effect of e ^q T operator on a two-body endshell HO state for 70hui included space. Plot a (c) shows the un-corrected 
HO state, while plot b (d) shows the ensuing endshell corrections. Note the large corrections at large-r. 



TABLE II: Triton B.E. at different oscillator 6s and As. These 
results ignore the small isospin-symmetry breaking effects in- 
cluded in the Avl8 potential. Hence isospin is conserved 
(T=l/2). Numbers in parentheses refer to the lowest eigen- 
value of the bare H within the corresponding model space. 
Results in MeV. 



A 


b=0.83 fm 


b=1.17 fm 


b=1.66 fm 





-7.6147 (218.533) 


-7.6144 (74.113) 


-7.6148 (26.480) 


2 


-7.6143 (65.260) 


-7.6144 (26.732) 


-7.6145 (12.556) 


4 


-7.6144 (26.226) 


-7.6141 (15.840) 


-7.6141 (9.068) 


6 


-7.6179 (13.268) 


-7.6140 (9.899) 


-7.6141 (6.524) 


8 


-7.6179 (5.985) 


-7.6136 (6.725) 


-7.6139 (5.257) 


10 


-7.6182 (1.606) 


-7.6131 (4.156) 


-7.6138 (4.093) 



TABLE III: Triton B.E. at different oscillator 6s and As. 
These results include isospin-symmetry breaking parts of the 
Avl8 potential. Hence the ground state consists of both 
T= 1/2 & 3/2 channels. Results in MeV. 



A 


b=0.83 fm 


b=1.17 fm 


b=1.66 fm 





-7.6219 


-7.6217 


-7.6213 


2 


-7.6230 


-7.6209 


-7.6209 


4 


-7.6255 


-7.6217 


-7.6206 





FIG. 2: Effect of e ^q T operator on a particular three-body endshell HO state for lOhus included space. Plots a and c show the uncorrected wavef unctions for the 1 So 
and 3 Si channels, respectively, while plots b and d show the ensuing corrections. Note the large corrections at small-p, corresponding to large corrections at large-r. 
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FIG. 3: Variational results for 3 He for various included spaces. 
The upper sets of lines correspond to calculations without 
endshell corrections, while the lower sets have endshell cor- 
rections. 




FIG. 4: WhuJ included-space calculation for 3 He using expan- 
sion given by Eq. ^3 



FIG. 5: Problematic configuration for three-body expansion. 
Such a configuration consists of two closely spaced particles 
(1 and 2), while a third particle is separated far from the first 
two. 
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FIG. 6: Diagrammatic expansion of e-qt-qv^ Q' ^ ne 
dashed lines correspond to V12 interactions. 




FIG. 7: Diagrammatic expansion of E _^q H Q f° r the three- 
body problem. The shaded ovals on the RHS are given by 
Fig. |S| Aside from the first term on the RHS, the remaining 



terms contribute to V, 



(3+0) 



I2,eff 



TABLE IV: Triton < 0hto\H e f f \0hu > matrix element at dif- 
ferent oscillator 6s and model sizes A. Within each cell, the 
top number represents < T e ff >, the middle number 

3 < ^/ 1 >j an d the bottom number 3 < V^t// >■ Num- 
bers are in MeV. 



hu A 





2 


4 


6 


8 


10 


00 


60 


64.51 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 




-68.36 


-38.41 


-32.54 


-26.40 


-21.03 


-16.52 


128.50 




-3.76 


8.26 


22.38 


37.25 


51.34 


64.06 


0.00 


30 


33.59 


45.00 


45.00 


45.00 


45.00 


45.00 


45.00 




-45.11 


-47.20 


-45.61 


-44.05 


-42.50 


-40.99 


29.10 




3.91 


3.71 


6.17 


8.96 


11.96 


15.03 


0.00 


15 


17.67 


22.50 


22.50 


22.50 


22.50 


22.50 


22.50 




-26.37 


-31.54 


-30.43 


-29.73 


-29.16 


-28.66 


3.98 




1.08 


2.43 


4.16 


5.03 


5.68 


6.26 


0.00 
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FIG. 8: Perturbative calculation of triton and He B.Es. for different As and bs using the expansion given in Eg. 1341 The solid 
black lines of each graph correspond to the exact binding energies. 
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FIG. 9: Relative sizes of T eff , SV^j], and 3^+°} for Qhuj 
and 2hw triton calculations as a function of 6. The upper pairs 
of lines refer to T e //, the middle pair V^^^and the bottom 

two Vjjf^Ji. Note that the respective sums of all three terms 
gives the triton binding energy. 
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TABLE V: A particular triton < 0huj\H e ff\2huj > matrix 
element at different oscillator bs and model sizes A. Within 
each cell, the top number represents 3 < V^.)\ >, and the 

bottom number 3 < V^^f >• Numbers are in MeV. Note 
that T e ff = for this case. 



hcu A 


2 


4 


6 


8 


10 


oo 


60 


-33.149 


-47.107 


-46.936 


-46.873 


-46.837 


-45.375 




0.146 


-1.771 


-1.692 


-1.016 


-0.330 


0.000 


30 


-13.400 


-18.898 


-18.772 


-18.665 


-18.564 


-18.952 




0.641 


-0.197 


-0.655 


-0.920 


-1.065 


0.000 


15 


-5.053 


-6.964 


-6.911 


-6.884 


-6.857 


-6.769 




0.398 


0.385 


0.296 


0.194 


0.104 


0.000 



TABLE VI: A particular Triton < 2hu)\H ef f \2huj > matrix 
element at different oscillator bs and model sizes A. Within 
each cell, the top number represents < T e // >, the middle 

number 3 < V^t// >, an <A the bottom number 3 < V^f t// >• 
Numbers are in MeV. 



hu> A 


2 


4 


6 


8 


10 


oo 


60 


122.210 


150.000 


150.000 


150.000 


150.000 


150.000 




-11.616 


11.714 


15.141 


18.365 


20.788 


89.884 




10.381 


18.639 


25.302 


31.312 


36.794 


0.000 


30 


61.986 


75.000 


75.000 


75.000 


75.000 


75.000 




-16.159 


-12.223 


-11.338 


-10.376 


-9.496 


21.034 




4.571 


7.211 


8.547 


9.790 


11.033 


0.000 


15 


31.679 


37.500 


37.500 


37.500 


37.500 


37.500 




-10.649 


-10.579 


-10.162 


-9.805 


-9.509 


4.019 




2.144 


3.365 


3.782 


4.071 


4.335 


0.000 



TABLE VII: Values of b where the contribution from V^J}) 
vanishes for a few included spaces sizes. Note that for each 
included space, there are two values. 



A 





2 


4 


b (fm) 


0.895, 1.838 


0.562, 2.219 


0.375, 2.706 



